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We present the second part of an if-band (1.6 fim) "atlas" of z < 0.3 SCR 
radio galaxies, using the Hubble Space Telescope Near Infrared Camera and Multi- 



\ Object Spectrometer (HST NICM0S2). We present new imaging for 21 recently 

Q\ \ acquired sources, and host galaxy modeling for the full sample of 101 (including 



11 archival) - an 87% completion rate. Two different modeling techniques are 
applied, following those adopted by the galaxy morphology and the quasar host 
galaxy communities. Results are compared, and found to be in excellent agree- 
^ \ ment, although the former breaks down in the case of strongly nucleated sources. 
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Companion sources are tabulated, and the presence of mergers, tidal features, 
dust disks and jets are catalogued. The tables form a catalogue for those inter- 
ested in the structural and morphological dust-free host galaxy properties of the 
SCR sample, and for comparison with morphological studies of quiescent galax- 
ies and quasar host galaxies. Host galaxy masses are estimated, and found to 
typically lie at around 2 x 10^^ Mq. In general, the population is found to be 
consistent with the local population of quiescent elliptical galaxies, but with a 
longer tail to low Sersic index, mainly consisting of low-redshift {z < 0.1) and 
low-radio-power (FR I) sources. A few unusually disky FR II host galaxies are 
picked out for further discussion. Nearby external sources are identified in the 
majority of our images, many of which we argue are likely to be companion galax- 
ies or merger remnants. The reduced NICMOS data are now publicly available 
from our website 0. 

Subject headings: galaxies: active — galaxies: fundamental parameters 



Introduction 



Extra-galactic radio sources have long posed an enigma to our understanding of the 
universe. While they can occupy some of the most gravitationally dominant galaxies in clus- 
ters, their bolometric luminosity is dominated not by the starlight from their galactic hosts 
but by the (often physically much vaster) regions of low-density radio-emitting material that 
have an origin in a tiny region of space at their very core. We now understand that the 
physics of these objects is related to the super-massive (> lO^M©) black-holes residing at 
their cores and the fueling of these black holes by accretion of gas and dust. However it 
is still unclear why so many apparently optically similar galaxies exist in the universe with 
no sign of radio activity. What m akes one black hole ra dio-active, and another not? The 
concept of a black- hole duty-cycle (IRichstone et al.lll998l ) as the fraction of black holes that 
are active at any time is a useful statistical characterization, but we still do not understand 
what physically drives it. It is known that, broadly, the comoving number density of black- 
holes evolves along with the cosmic star-formation rate, and this can be seen as a correlation 
betw een the evolution of the black-hole duty cycle, and that of the cosmic star formation 
rate (IWang et al.ll2006l ). Mergers have long been regarded as a strong candidate for trig- 
gering AGN activity, providing a mechanism for the delivering th e necessary (gas and dust) 
fuel to the central regions of a galaxy and the black-hole therein (IKauffmann and Haehnelt 



^http: / /archive. stsci.edu/prepds/3cr/ 
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2000l : iDi Matteo et al.l l2005l ) . Studies have shown the presence of nuclear dust disks i n the 



cores of many radio galaxies (ISadler and Gerhard! Il985 



van 



Dokkum and Franx Il995 ). but 



similar numbers are also found in quiescent galaxies (IVeron-Cetty and Verorull988l ). Early 
studies of quasars showed morphological dis t urbances or tidal features t o be present in mari y 
of the host galaxies (e.g. ISmith et al.lll986l : iHutchings and Nefllll992l : iBahcall et al.l 119971 ). 
but more recent (and better-selected) samples have shown that the proportion exhibiting 
majo r tidal features is indistinguishable fro m that in the quiescent elliptical galaxy popula- 
tion (IDunlop et al.ll2003l : iFloyd et al.l l2004j ) . Thus today it remains unclear whether there 
are any optical-IR properties of radio galaxies that distinguish them as a class from the pu- 
tative parent population of quiescent elliptical galaxies, rather than simply being inherited 
from them. 

In recent years, much attention has been drawn to the correlatio n seen locally between 



bulge and black hole mass i n normal quiescent elliptical galaxies (e.g. iGebhardt et al.ll2000 



Merritt and Ferraresell200ll . AGN provide two obvious possible feedback mechanisms under 
which such a relation might arise, through radiation pressure and through the action of ra- 
dio jets emanating from the central eng ine. Theoretical de velopments, and recent numerical 
studies using semi-analytic models (e.g. ICroton et al.ll2006l ) to explore the evolution of galax- 
ies have shown that AGN feedback can in principal enforce such a correlation by shutting off 
star-formation and ejecting gas and dust from the most massive galaxies. The same models 
offer a natural solution to the problems of over-production of very massi ve galaxy haloes, and 



simultaneously can explain the so-called "re d sequence" o f galaxies (e.g. ICattaneo et al.ll2006 



Pi Matteo et al.ll2005l : ISpringel et all 120051 ). Best et al.l fl2006r ) showed that the radio-loud 
phase of AGN activity is able to suppress the cooling of the host galaxy halo gas sufficiently 
that the radio source can control the rate of growth of the elliptical galaxy. 

It is important to investigate whether any statistical differences can be detected between 
samples of radio-loud AGN and samples of quiescent galaxies or radio-quiet AGN. One 
difficulty is the way in which active and quiescent galaxies have traditionally been studied 
is rather different. The most detailed morphological studies of large samples of quiescent 
galaxies are not conducted in the same way as studies of radio galaxies, or the host galaxies 
of quasars and other AGN. To investigate the causes of AGN activity and radio-loudness 
we need large samples of objects in which both AGN and environment can be explored at 
multiple wavelengths. The 3C makes an ideal such sample for studying the effect of radio- 
loud AGN activity on host galaxy environment, and vice- versa. Here we focus on t he 3CR, 



a wel l -defined sub-sample of the most powerful northern hemisphere radio galaxies (IBennett 



19621 : ISpinrad et al.lll985l ). We aim to provide a complete inventory of the properties of the 
low-redshift {z < 0.3) 3CR, including physical characterization of the host galaxies free of 
dust to determine the underlying galaxy structure and dominant stellar mass; the presence of 
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companions, mergers, dust disks, and jets; and the clear detection of any unresolved optical- 
IR nuclear sources. This will enable us to explore whether anything distinguishes them from 
normal quiescent galaxies, and continue the ongoing study into the question of what makes 
radio-galaxies "radio loud". 



1.1. The SCR snapshot program 

In HST cycle 13 (2004-5) we embarked on a Near-Infrared Camera and Multi-Object 
Spectrometer (NIC MOS) if-band snapshot imaging campaign of z < 0.3 3CR sources 



( iMadrid et al.ll2006l - hereafter paper I) . The infrared is an important realm for our grow- 
ing understanding of the relationship between AGN and galaxies, enabling us to study the 
underlying host galaxies free from the distorting effects of dust. The SNAP program is now 
complete, having continued through cycle 14 (2005-6) at reduced priority. Since publica- 
tion of paper I (which presented images for the first 69 targets, observed during cycle 13) 
an additional 22 sources have been successfully observed, and images for these new targets 
are presented here. Finally, archival data for an additional 11 objects previously observed 
with NICM0S2 in F160W were obtained from the Multi-mission Archive at Space Telescope 
(MASli). 

In this paper we seek to characterize the local 3CR host galaxy population in such a way 
that they can easily be compared to similarly large samples of quiescent galaxies, mergers, 
radio-loud and radio-quiet quasar host galaxies, and brightest cluster members. We present 
modeling of all the data using elliptical isophote (1-dimensional radial profile fitting) and two- 
component 2-dimensional galaxy-modeling techniques. We present the sample, observations 
and data reduction in section [21 We discuss the ID and 2D modeling of the NIR {H- 
band) host galaxies in section [31 In section [H we present notes on, and images of, the newly 
observed and archival targets that were not presented in paper I. Basic host galaxy properties 
are presented, and the modeling techniques compared in section [51 In section [61 we present a 
census of the companion sources and merger environments of the sample. Section [7] presents 
a general discussion of the sample and study and suggestions for the future. Section [HI 
concludes with a summary of our main findings. Detailed comparison of the sample with 
control samples of ellipticals, quasars and mergers is left to a companion paper (Floyd et al. 
in preparation). We have made all of our reduced data publicly available on the internet§|. 



^http: / / archive.stsci.edu/ 

■^http: / / archive.stsci.edu/prepds/3cr/ 
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2. Sample and data reduction 

The majority of tlie data tliat we use in this paper were taken during the HST Snapshot 
Program SNAP-10173 (PI: Sparks). The near-infrared images, fluxes and notes for the first 
69 of our targets were presented in paper I. In this paper we present the imaging for 21 
additional targets that were unobserved at the time of paper I's publication. An observing 
log for these newly observed sources is presented in table [H In addition, we analyze archival 
NICM0S2 F160W data that exist for a further 11 objects - see table H Note that NICl and 
NIC3 lack the field of view and the resolution (respectively) for this morphological study, 
so observations of SCR sources on those two chips in the archive have been omitted from 
study here. Finally, the NIC2 observations of 3C 273 were unsuitable for the present study 
as they provide only 30 s of on-galaxy integration, with the bulk of the observing time being 
dedicated to the famous jet. Altogether this selection produces a near-complete sample of 
101 powerful {L^ghz < 10~^^ W Hz~^ sr~^) northern hemisphere radio sources at z < 0.3 
observed in a single infrared band. Radio luminosities (at 178 MHz) for the full sample are 
plotted in Fig. [TJ Other basic sample data are presented in table [31 Notes on, and images of 
the newly observed and archival targets are presented in section HI See paper I for images 
of, and notes on the other 69 sources. 

We have re-reduced all of the original data presented in paper I, along with the new 
observations presented here, as described below. This new technique conserves the sky back- 
ground throughout for all sources (allowing for better determination of the image statistics 
which is essential for accurate morphological characterization), and maintains the original 
pixel scale of 0'.'075, which is near-critically resolved, and perfectly adequate for large-scale 
morphological studies such as this. On a large sample, the modeling of dithered images with 
double the spatial resolution results in an impractical increase in the computational time 
due to the convolution of each image with a point spread function. 

We retrieved the data, pre-processed through the standard calibration pipeline, from the 
MAST. We performed the data reduction described below using the Space Telescope Data 
Analysis System (STSDAS) software running under PyRAF. Two anomalies not removed by 
the calibration pipehne are corrected during our data reduction: the pedestal effe ct, and the 



coronographic hole. NIC2 has a quadrant offset bias known as the pedestal effect (INoU et al. 



20041 ). We use the task pedsub to remove this source of error, pedsub eliminates the 
pedestal effect but leaves the sky level untouched. We found that using pedsky on nearby 
sources (as in paper I) led to some errors in the photometry due to inaccurate characterization 
and subtraction of the background flux. The NICMOS coronographic spot appears as a small, 
well-defined region of erroneous flux on each individual dithered image, and is simply masked 
out before combining the images. 
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Each object was observed using a square dither pattern with four exposures of 288 s 
each (1152 s total integration time on each source). We use Multi drizzle, a one-step tas k 



to perform image combination, dithering and cosmic ray rejection (iKoekemoer et aLll2002l ). 
We set Multidrizzle to leave the NIC2 native pixel size of C/OTS and to not perform sky 
subtraction. Note that Multidrizzle produces a final output image in units of electrons 
per second. We convert to DN before proceeding, by dividing the images by the Analogue-to- 
Digital Conversion (ADC) gain and multiplying by the integration time, obtained from the 
relevant FITS image headers (ADCGAIN and EXPTIME respectively). Also note that the 
default behavior of multidrizzle is to output weight maps that are exposure time maps. 
We adjusted this by setting final_wht_type = "ERR" to produce inverse-variance maps. 

For each source we have fc-corrected the apparent magnitudes to rest-frame H-band as- 
suming a spectral index of a = 1.5 for the host galaxy and a = 0.2 for the nucl eus where nec 



essary {fy oc z/ "). Galactic extinction corrections have been applied following ISchlegel et al. 



(Il998l ). We assume throughout a flat, A-dominated cosmology in which Q\ = 0.7, and 



ho = 0.7. 



3. Modeling Procedure 



We have used Ellipse ( jJedrzejewskil 119871 ) and Galfit (jPeng et al.ll2002l ) to model each 
source in the sample. The steps required to implement these models properly are outlined 
below (section [3.ip . To the elliptical isophotal values returned by Ellipse, we fit 1-dimensional 
Sersic -law quasi-radial profiles (section [3l2!) . Properties of these profiles are compared to 
those obtained by direct 2-dimensional fitting to the image data usin g Galfit (section 13.31) . 
We make use of the summary of photometric galaxy relationships by Caon et al. ( 19931 ) to 
convert between surface brightness, scale-length and total flux for different values of Sersic 
index, n. 



3.1. Pre-modeling 



3.1.1. Masks 



We began by creating a mask and error frame for each target, using the output weight 
frame and hot pixel flags from MULTIDRIZZLE. The input bad pixel masks are combined 
by Multidrizzle to produce a final mask of all hot or bad pixels. To this we added an 
additional source mask to remove any external sources, companions, diffraction spikes, and 
the warm (ampglow- affected) corners of the image from the fit. These masks were applied 
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Fig. 1. — 178 MHz radio luminosity of the sample. The dotted line indicates the conventional 
radio- loud / radio-quiet cut-off for quasars (L5 qhz = 10"^'' W Hz~^ sr ), converted to 
178 MHz assuming a power-law index of 1 for the radio luminosity (Lj, ~ ^^^)- The solid 
line indicates the (9 Jy at 178 MHz) flux limit of the SCR survey. Note that two of the most 
well-known low-redshift targets, 3C 71 (M 77, NGC 1068), and 3C 272.1 (M 84, NGC 4374) 
actually fall below the radio-loud quasar cut-off. Upward-pointing triangles represent FR II's; 
Downward-pointing triangles are FR I's; Diamonds are High Excitation Galaxies (HEG's); 
Squares Low-Excitation Galaxies (LEG's). The single solid star represents 3C 71. 
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to all the subsequent fits described below. 



3.1.2. Errors 

The combined weight (l/c^) frames that are output by Multidrizzle incorporate the 
Poisson noise and read noise on each pixel. We converted them to standard error (cr) frames. 
For the nuclear-dominated sources (with nuclei providing > 50% of the total flux), the errors 
are dominated not by Poisson or read noise, but by the "sampling error" associated with the 
lack of know ledge of the point spread function (PSF). In these cases we follow the technique 



described in lFloyd et al.l (120041 ). to provide accurate representation and weighting across the 
entire dynamic range, from the nucleus to the wings of the PSF. In practice, this means 
estimating the "sampling error" in concentric circular annuli, centered on the nuclei, from 
the deviation of the source from the PSF within each annulus. Where this drops to the 
Poissonian level (typically within a radius of ~ O'.'S), we revert to the properly calculated 
individual pixel errors described above. 



3.1.3. Point spread functions 



Point spread functions (PSF's) were calculated for each target using Tinytim (IKrist 



19991 ) ■ setting the position of the PSF to the average of the four central positions of the target. 



We generated a PSF at twice the resolution of the NIC2 chip, to provide oversampling, and 
accurate shifting and convolution. Defocus was estimated for the time of each observation 
using the tables and figures provided on the NICMOS website 0, and this information was 
included in each Tinytim model. The Tinytim models were, last of all, rotated through 
the same angle as the target images, to give north up in each case. 



3.1.4. Sky background level 

For all targets, an initial measurement of the approximate background level in the image 
was made using imstatistics to determine the photometric statistics on several empty 
regions of the chip. We attempted to measure the background in all four quadrants of the 
chip, as a check of the success of the quadrant bias removal. In practice, on all images 
where this was possible, the quadrants were found to have identical backgrounds, within 



"^http: / /www. stsci.edu/hst/observatory/focus/documents. html 
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the Poissonian errors. On objects that are too close for us to observe any background on 
the small 19 x 19" NIC2 chip (in practice, almost all objects at z < 0.1 and some large 
object s at higher redshifts) we employed the technique developed by ISparks and Jorgensen 
(119931 ): All of these targets have been observed previously using WFPC2 in either R or V 
band (or both). We obtained the reduced full WFPC2 mosaics from the SCR database 0, 
and measured the background flux on these larger images. Profiles of the flux in each of 
the 2 bands were produced, and plotted against each other. In all cases there is a region of 
linearity in the 2 flux levels, and a linear fit was made in this region, extrapolating down 
to obtain the IR flux at the level of the optical (WFPC2) background. As a check of the 
robustness of this fitting technique, we also deployed it on images of five higher redshift 
sources, where a clear measurement of the infrared background was also possible from the 
NIC2 image. In these cases we produced a small cutout of the infrared image, such that no 
background light was visible, mimicking the situation for the lower redshift targets. In these 
tests, all five fits yielded infrared backgrounds that, again were within the Poissonian error 
of the true, measured background fluxes on the original infrared images. 



3.2. Ellipse fitting 

We used the Ellipse modeling procedure in Iraf.stsdas to fit elliptical isophotes to 
the sky-subtracted NICMOS images. We initialized each fit by feeding it the position of 
the centroid of the source, as measured using imcntr, and estimated values of the position 
angle (PA) and ellipticity (1 — b/a). We began the fit at 10 pixels (semi-major axis length) 
working outwards in geometric steps, increasing in size by a factor xl.l. Once the fit has 
terminated due to insufficient pixels (less than half on a given annulus, fflag = 0.5), or 
convergence (largest harmonic amplitude < 5% of the residual RMS around the harmonic 
fit) it returns to fit the inner region with identical geometric steps, all the way to the centre 
of the target. The central position, ellipticity and PA of each ellipse were left free to explore 
the stability of the isophotal centering, and investigate structures in the galaxy deviating 
from the purely elliptical. Deviation from a perfect ellipse at each radius is broken down 
with a Fourier series, characterized by the 3rd and 4th order sine and cosine coefficients. 
We convert from the E LLIPSE parameters A3, A4, B3 and B 4 to 04/0, 64/6 following the 



conversion described in iMilvang- Jensen and J0rgenseru (119991 ) - see that paper for a useful 



discussion of all the ellipse fit parameters. Additional masking was performed, in software, 
on pixels deviating by more than Sa^nnuius from the mean value on a given fitted elliptical 
annulus. 



''^http://acs.pha.jhu.cdu/^niartel/ 
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Modeling with ellipse allows us to track radial changes in the ellipticity and position 
angle of the isophotes of galaxies, providing a powerful analysis of the structure of the galaxy. 



and p resence of any disturbances from a relaxed elliptical morphology. 1-dimensional iSersic 
(119681 ) profile fits were made to the best-fitting ellipse-isophotal intensities, after convolution 
of the model with a point spread func tion (PSF), using a Levenberg-Marquardt non- linear 
least-squares minimization technique (jPress et al.lll992l ). This techniq ue has been widely 



used in the past fo r studying the pr operties of elliptical galaxie s (e.g. iBender et al.l Il992 



Graham et al.ll 19961 ) and mergers (e.g. iRothberg and Josephll2004l ). but note that it does not 



provide a true "radial profile" since adjacent elliptical isophotes may have different position 
angles and ellipticities, thus overlapping, and contributing flux to one-another. For purposes 
of comparison with the 2D fits (see below) and with the literature, the median values of 
ellipticity, position angle, and diskiness were computed for the isophotal fits between a radius 
of O'.'S (inside of which the nucleus may dominate) and where the signal drops to Sa^ky, or 
the edge of the chip if that is reached first. The results of the elliptical isophote fitting are 
presented in table HI 



3.3. Galfit 



Galfit (jPeng et al.l |2002| ) is a versatile fitting code which allows the user to fit multi- 
ple components to a galaxy in an iterative fashion, using a downhill gradient Levenberg- 
Marquardt technique. The model galaxy is convolved with a PSF at each iteration, and 
compared to the data through the statistic. We fitted each object using a single Sersic 
component convolved with a PSF, plus an additional unresolved nuclear component rep- 
resented by the same PSF. Error frames and masks were produced as described in sec- 
tion |^^2_§ibove. The ma jority of objects are too distant to allow detection of a "Nuker" 



[e.g. iFaber et al.l 119971 ) or I Graham et al.l (120031 ) type core, so we adopt just a single Sersic 
component for study throughout the sample for consistency, simplicity, and comparison with 
the bulk of the literature. 

The complete (Sersic -|- PSF + sky) model has 11 free parameters: Centre of host {x,y); 
Luminosity of the host galaxy (Lh)', Scale length of the host galaxy (Re); Position Angle 
of the host galaxy (9); Ellipticity of the host galaxy (b/a); Sersic parameter of the host 
galaxy (n); Luminosity of the nucleus {L^); Position of nucleus {x,y); Background flux. We 
initially held fixed the values of the Sersic index (n = 4 - i.e. a de Vaucouleurs elliptical), 
the "diskiness" parameter, (c = 0), and the sky background flux (at a level determined 
as described in section [3.1.41 above). After a satisfactory fit was obtained, we freed these 
parameters one-by-one, in an iterative modeling procedure, in the order listed above. Last 
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of all, the sky value was freed up to see if a better fit could be obtained, and to test the 
degeneracy between sky background and host galaxy flux. The models presented here include 
this fitted sky value, which was always consistent with the measured sky value. The results 
of the galfit Sersic model fitting are presented in table O 



As discussed above (section [3X2]), we have adopted the error analysis technique of lFloyd et al 



(120041 ) for nuclear-dominated objects. Since this is not the standard error treatment used 
by Galfit, we performed extensive tests (discussed in the Appendix) of the effect of this 
assumption on the fits produced by Galfit using both real and synthetic sources. We found 
Galfit to prod uce identical result s using each technique for sources with weaker nuclei. We 
argue that the iFloyd et al.l (120041 ) technique is the formally correct analysis of the errors for 
sources with dominant nuclei since it properly accounts for all sources of error across the 
image. 



3.4. Errors on the fitted parameters 

Errors on the fitted parameters in each case are calculated from the diagonal elements of 
covariance matrix for that set of parameters, following the assumption that each parameter 
is independent. The errors quoted in tables and H] are formal l-cx statistical errors. Note 
that in many cases, the true errors are dominated by the systematics, since the galaxies are 
not ideal. To investigate this, we performed simple tests of different galaxy fitting algorithms 
as described in the appendix. Including all systematics, we expect that all our morphological 
results are accurate to within 10%, and host fluxes within 2% 



4. Notes on individual objects 

In this section we present notes and images for the newly observed objects in our SNAP 
program (i.e. those not included in paper I) - Fig. [21 and those taken from the archive - 
Fig. [31 



4.1. Newly observed sources 

4.1.1. 3C 15; z = 0.073 



This source exhibits a prominent optical-IR jet (iMartel et al.lll998l ) in an apparently 
undisturbed round elliptical that almost fills the NICM0S2 field-of-view. No bright nuclear 



point source is detectable in H-band. 



4.1.2. 3C 17; z = 0.219 



This very peculiar radio source (iMorganti et al.l Il999l ) and BLRG has a quite round 
elliptical galaxy, and a prominent IR nucleus. There are several faint companions in the 
immediate surroundings of the host galaxy. 



4.1.3. 3C 33; z = 0.059 

An FR II HEG source with a prominent dust disk aligned ~ 45° to the radio axis. The 
host galaxy is elliptical, appearing slightly disturbed due to the dust disk. 



4.1.4. 3C 98; z=0.030 



A NLRG FR I, hosted by an undisturbed elliptical with smaller companion elliptical. 
No internal structure is visi ble in the NI CM0S2 image. The radio source is double lobed, 
running almost north-south (jMillerlll985l ). 



4.1.5. 3C 132; z = 0.2U 

3C 132 appears near the western corner of the NIC2 image, and has a somewhat elon- 
gated, smooth elliptical host galaxy. Several companions are visible, including one bright 
foreground star at the image center. The radio axis of this FR II source is almost perpen- 
dicular to the major axis of the host galaxy. 



4.1.6. 3C 153; z = 0.277 



3C 153 is a radio galaxy (angular size 6 arcsec=25 kpc) oriented at P A^ 50° with 
an FRII radio morphology and a large arm-length ratio ~ 1.9 (ILaingl Il98ll ). It has been 
classified as a CSS source (e.g. lAkujor et al.l Il99ll ) but it would be at the upper end of 
the size distribution for such sources. The [ Oil] emission l ine n ebula is slightly extended 
along the radio axis (IMcCarthy et al.lll995l ). Ide Koff et al.l (Il996l ) note that the galaxy has 
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an elliptical nucleus and strong emission lines. The NICMOS image shows a fairly round 
host galaxy. There are a couple of unresolved sources nearby but they don't correlate with 
knots or hotspots in the radio source. Two other features are caused by the intersection of 
diffraction spikes. 



J^.l.l. 3C 166; z = 0.245 



An unusual radio source featuring two lobes with very different morphologies (ISpangler and Bridle 



19821 ). The radio axis runs north-south, with Spangler & Bridle's southern features C and D 
detected on the IR image as infrared hotspots. The source has a clearly-detected, unresolved 
infrared nucleus. 



4.1.8. 3C 234; z = 0.185 



3C 234 exhibits a brig ht quasar-like nucleus in the NICMOS image, and is known to be 
a broad-line ra dio galaxv (jMcCarthv et al.lll995l) with classical double FR II radio source. 



In the optical (jde Koff et al 



19961: 



McLure et al.lll999l ) there are features emanating from 



the nucleus in the east and west direction. The western feature, described as a tidal arm by 
McLure et al. is detectable on the new NICMOS image. There are two additional sources to 
the west in the same direction, the furthest of which is also detectable in F555W, and F702W. 
Both are detectable in F675W. This galaxy is known to have strong extended emission lines 
which likely explain these features. Pa f3 (1.28 /xm) is in the passband of F160W at the 
redshift of 3C 234. 



4.1.9. 3C 258; z = 0.1Q5 



This galaxy is a well-known compact steep spectrum (CSS) radio source, with a com- 
pac t, double radi o struc t ure seen on scales of a few tenths of a n arcsecond along PA:^ 



33° ISpencer et al.l (119891 ): lAkujor et al.l (Il99ll ): iLudke et al.l (119981). although its arcminute 
scale structure is aligned along a i iearly north- s outh direction IStrom et al.l (Il990l ). Its red- 
shift h as been quoted as z=0.165 ISmith et al.l (119761 ) . but this value is in dispute as iDey 
(119941 ) note that spectra taken at Lick show a likely background quasar, based on the obser- 
vation of a single emission line at 711lAwhich would give z = 1.54 if due to Mg II 2798A. 
However, due to poor seeing at the time the spectra were taken, no firm claims can be made 
as to the origin of this line. Our NICMOS data are consistent with a significantly higher 
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redshift, as they show a very compact host galaxy that is much fainter than one would expect 
for an object at z=0.165. The nuclear point source in our image is also quite bright, with an 
arc-like object 3" ESE of the AGN which may be an irregular foreground {z = 0.165) galaxy 
causing the confusion. 



4.1.10. 3C 284; z = 0.239 

An FR II HEG radio source, hosted by a disturbed elliptical with a SE tidal tail toward 
its most prominent companion 4"SE. Several fainter companions are also visible on the 
NICMOS chip. In the optical the source has a disky appearance and double nucleus, owing 
to this same tidal tail and associated dust lanes. 



4.1.11. 3C 296; z = 0.025 



3C 296 is hosted by an extremely boxy elliptical galaxy. It h as a double lobed j et with 
clear radio emission along the jet axis (FR I). The WFPC2 image ( iMartel et al.lll999l ) shows 
a truncated edge-on disk embedded in the galaxy, similar to that seen in NGC 4261. The 
disk is aligned with the major axis of the elliptical galaxy, with the jet perpendicular to the 
disk. The nucleus is seen in both WFPC2 image (though faint) and in the NICMOS image 
(brighter). X-ray and radio observations of this galaxy are discussed in iHardcastle et al. 
(120051 ). 



4.1.12. 3C 300; z = 0.270 

A double-sided FR II HEG, hosted by a compact, elongated ~ 3 kpc elliptical host 
galaxy with a faint companion 8"East, and a very faint tidal distortion. The radio axis runs 
roughly SE-NW, and there is a faint thread of emission along this axis on the NW side, 
which is a candidate jet. However, it is almost lost in the noise of the NICMOS chip edge. 



4.1.13. 3C 323.1; z = 0.264 

This QSO has a bright nucleus that produces diffraction spikes in the near-infrared 
image. It has a close companion to the northwest that does not show clear signs of in- 
teraction, but is consistent with being the remnant nucleus of a galaxy in the final stages 
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of merging (IMiller and Sheinia l2003l : ICanalizo and StocktonI 119971 ) . The Keck spectroscopy 
presented by Miller et al. demonstrate that the host is not a normal elliptical galaxy, and 
from our if-band image there is some evidence of disturbance, which may explain the small 
scale length fit to the data. It is also possible that the radio emission from this source is 
Doppler boosted, making it an artificially radio-loud object. 



4.1.14. 3C 326; z = 0.090 



A double source - t he new NICMOS imag e shows only the southern, which is the stronger 
of the two in the radio (jRawlings et al. 1991). A larger companion elliptical sits 25" (43 kpc) 
to the north. A dust disk bisects the nucleus of the small elliptical host galaxy, giving the 
appearance of a disk. 3 very faint small companion sources are visible within, or close to the 
main galaxy halo. Unusually the radio axis is aligned close to the axis of the disk. 



4.1.15. 3C 332; z = 0.270 

A powerful FR II radio galaxy with a prominent quasar-like infrared nucleus marking 
the center of an elliptical host galaxy. There is a smaller companion 3"to the SW. 



4. 1.16. 3C357;z = 0.im 



3 C 357 is ail FR I I radio galaxy with angular si ze 90"(250 kpc) at a PA 110 de- 
grees ( Fanti et al. 1987 : Harvanek and Hardcastlj Il998l ). In both lobes the hotspots are 
off to one side suggesting there has been a change in the jet directioi i. There is emission liri e 
gas extending a few arcse c from the nu c leus al ong the radio axis (IMcCarthy et al.l Il995l ) . 
de Koff et al.l (119961 ) and ICapetti et al.l (l2000al ) note that there are filamentary dust lanes 
southwest of the nucleus. The galaxy isophotes are very elliptical with no remarkable fea- 
tures. There is a single faint candidate companion source 8"East, with a larger companion 
just on the eastern edge of the chip. No unresolved nucleus is detected, but the core of the 
galaxy appears unusually bright, and the host has a high Sersic index {n = 5.34). 
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4.1.17. 3C 403.1; z = 0.055 

This double radio source lies at low galactic latitude {b ~ —14°), resulting in numerous 
stars on the image. In spite of this one can see that the undisturbed-looking elliptical 
host galaxy is located in a fairly dense environment, with one large com panion seen on the 



NICMOS chip, and 3 others similar in brightness visible on the WFPC2 (IMartel et al.lll999l ) 
image. 



4.1.18. 3C 410; z = 0.248 

Our observations of 3C 410 suffered a pointing error, resulting in only two usable images 
of the source - hence the rather noisy image. A strong central nuclear point source is visible 
centered on a round elliptical host galaxy. Numerous unresolved sources are visible in the 
field. 



4.1.19. 3C 424; z = 0.127 
Faint tidal distortions are visible to the north on our NICMOS image, and t his source lies 



in a d ense environment, with numerous large companions visible on the WFPC2 (jde Koff et al 



19961 ) image. The radio source apparently lies on the edge of the cluster. Two bright un- 
resolved sources are visible to the north, with a resolved object NW, on the edge of the 
NICMOS chip, barely detected at R. 



4.1.20. 3C 442; z = 0.02Q 



Only two exposures can be used of the four, due to telescope pointing err ors. The host 



galax y is a smooth elliptical, with 2 companions visible on the WFPC2 image (IMartel et al. 



19991 ). The point source on the edge of the NICMOS chip to the WSW is in line with the 



radio axis, and is a candidate hotspot - also clearly visible on the WFPC2 images. 



4.1.21. 3C 459; z = 0.219 



A strong IR nucleus, and heavily distorted off-center host galax y distinguish this FR I I 
radio galaxy. The host is dominated by young stellar population (ITadhunter et al.ll2002l ). 
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and must ha ve undergone a recen t merger. The radio source is small (10") and also very 
asymmetric (IMorganti et al.lll999l ). 



- 18 - 




Fig. 2. — NICM0S2 fl60w [H band) images of the newly observed objects aligned north up 
- see table [H 
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Fig. 2. — (continued). NICM0S2 fl60w {H band) images of the newly observed objects 
aligned north up - see table [H 
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Fig. 2. — (continued). NICM0S2 fl60w [H band) images of the newly observed objects 
aligned north up - see table [H 
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Fig. 2. — (continued). NICM0S2 fl60w (H band) images of the newly observed objects 
aligned north up - see table [H 
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4.2. Archival data 



4.2.1. 3C 71 (NGC 1068, M 77); z = 0.003793 



A well-known spiral galaxy, and the archetypal Seyfert 2, NGC 10 68 is not truly "radio- 
loud" in terms of its radio luminosity. It has an obscured Sy 1 spectrum ( lAntonucci and Miller 
19851 ) and a compact (few arcsecond) radio jet a ligned NNE-SSW. The NICM OS data were 
obtained as part of the NICMOS GTO time by [Thompson and CorbinI (119991 ). The galaxy 
dwarfs the small NIC2 chip, and even the WFPC2 mosaic (PROPOSID's 5479, 5754). We 
were unable to produce an accurate photometric fit to this object, since we could not esti- 
mate the background flux. The spiral structure and active nucleus are readily apparent on 
the NICMOS image, as are numerous globular clusters. There is a vast literature on this 



source. 



4.2.2. 3C 84 (NGC 1275, Per A); z = 0.017559 



The central galaxy of the Perseus cluster contains an FR I radio source and Seyfert 2 
nucleus at the centre of a cooling flow. It shows a large number of companion sources in the 
infrared, and a complicated dust ni orphology in the optical-IR. The NICMOS data was first 
published by iMartini et al.l (120031). Only the centr al region is visible in the NIC2 image. 
We used the F702W image from iMartel et al.l (119991 ) to estimate the background. However, 
due to the immense size of this cD galaxy, even the full WFPC2 mosaic (PROPOSID 6228) 
does not detect the edge of the galaxy, and the results are photometrically uncertain. 



4.2.3. 3C 264; z = 0.021718 



3C 264 is a tailed radio source in the cluster Abe ll 1367§. On the a r csec scale there i s 
a bright one-sided jet wh ich is detect ed in the radio (jSaum et al. 



optical (ICrane et al.lll993l ). and X-ray (jTilakll2006l : iPadgett et al 



200 



19971 

T 



Lara et al. 



Baum et al 



1999), 



h997f ) 

note the existence of an inner dusty disk which causes an apparent ring in the galaxy surface 
brightness where the disk ends. The NICMOS image shows a brig ht nucleus, the jet and the 
ring noted by Baum et al. The data was previously published by lCapetti et al.l (l2000bl ). 



http://www.jb.man.ac.uk/atlas/dragns.htnil 
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4.2.4- 3C 270 (NGC 4261); z = 0.007465 
A well-studied nearby F R I with a boxy elliptical host galaxy and a circumnuclear dust 



disk. iFerrarese et al.l (Il996l ) used imaging and spectroscopy of this source to determine its 



black-hole m ass. A faint unreso lved nuclear point source is visible in the IR. Data first 
published bv lOuillen eTaP J200oh . 



4.2.5. 3C 212.1 (NGC 4374, M 84); z = 0.003536 

A large well-studied elliptical in the Virgo cluster with a prominen t dust lane visible 
in the optical, and faintly detectable in the NICMOS image iBower et al.l (120001 ). The radio 
source is a double-sided FR I jet and core. However, its radio power is somewhat lower 
th an the fiducial cu t off fo r radio-loud quasars. It is described as a "weak radio galaxy" 
by lLaing and Bridld (119871 ). 



4.2.6. 3C 274 (NGC 4486, M 87, Virgo A); z = 0.004360 

An FR I with ce lebrated optical jet. The NICMOS data was originally presented 
by (jCorbin et al.ll2002l ). where the source is studied in detail. 



4.2.7. 3C 293 (UGC 08782); z = 0.045034 

A merger-remnant with significant distortion at optical and IR wavelengths. The promi- 
nent optical dust lanes partiall y cover an op t ical jet that is clearly seen in the infrared. This 
object was studied in detail by lFloyd et al.l (l2006bl ). 



4.2.8. 3C 305 (IC 1065; UGC 09553); z = 0.041639 

A nearby FR I and merg er remnant, with prom inent dust lane across the nucleus in the 
optical. Studied in detail by IJackson et al.l (120031 ). The host galaxy is a large, elongated 
elliptical, with clear signs of disturbance in the optical and the IR. See lJackson et al.l (120031 ) 
for original publication of the data, and a detailed study of the object. 
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4-2.9. 3C 317 (Abell 2052; UGC 09799); z = 0.034457 



3C 317 is a slightly elongated elliptical with an unresolved nucleus. Two small elliptical 
galaxies are present in the field of view, one is located 7" to t he northeast, and th e second 
one 9" to the northwest. The image of this galaxy taken by iMartel et al.l ( 120021 ) shows a 
peculiar UV filament ~ 4 kpc south of the nucleus. This fila ment is described as a region of 
active star formation likely triggered by a recent merger by iMartel et al.l (120021 ). The data 
was first published by iQuillen et al.l (120001 ) . 



4.2.10. 3C 338 (Abell 2199; NGC 6166); z = 0.030354 



An FR I radio source in a dense environment. Two companions are prominent on the 
NICMOS image. The elliptical host galaxy is disturbed, with a dust lane crossing the nucleus 



faintly visible in the IR. First published by I Jensen et al.l (120011 ). See also lRavindranath et al. 

kooih . 



4.2.11. 3C 405 (Cyg A); z = 0.056075 



3C 405 is the archetypal FR II r adio galaxv (|Carilli and Barthellll996l ). There is a buried 
quasar detecte d in the host galax y ( Djorgovski et al.[l991 ). The hot spots are detected in 
the X-ray (e.g. IWilson et al.ll2000l . The galaxy envelope fills the 19 x 19"NIC2 chip, with the 
hot spots lying outside the NICMOS image. The inner galaxy contains interesting patchy 
and filamentary structure. There is strong dust absorption running roughly E-W below the 
nucleus. There is a strong patch of dust just NE of the nucleus. There is also what looks 
like part of an ionization cone pointing to the north-west. At the highest brightness levels, 
we see a bright point source centered on a rough "X" shap e which may define the edges of 
the ionization cone. Similar results have been described by iTadhunter et al.l (119991 ). 
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Fig. 3.— NICM0S2 fl60w 



[H band) images of the archival objects - see table [2j 
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Fig. 3. — (continued). NICM0S2 fl60w {H band) images of the archival objects - see 
table H 
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5. Host galaxy modeling results 

In this section we present results for the 1-D and 2-D fits, a basic comparison of the 
results from the two techniques, and statistics describing the modelled properties of the host 
galaxies, for comparison to other galaxy samples. 

5.1. Ellipse quasi-radial profile fits 

Properties of the ellipse fits are presented in table HI We note that the ellipse-fitting 
technique is only reliable in cases where there is no strong nuclear point source. In the 
presence of a strong nuclear point source, much of the galaxy structure is obscured by the 
asymmetric PSF, the central region of the galaxy is entirely obscured, and thus the galaxy 
cannot be accurately fit. In these cases, the 2-dimensional fit is able to perform far better 
as it can take account of the asymmetrical 2-D structure of the PSF. We also note once 
again, that the 1-D approach does not provide a true "radial profile" since adjacent elliptical 
isophotes will in general have different position angles and ellipticities, thus overlapping and 
contributing flux to each other. 

5.2. Parametric galaxy fits 

The direct 2-D Sersic fits are summarized in table [51 The technique is unable to follow 
all of the isophotal twists in a disturbed galaxy's profile, but provides an effective tool for 
modeling the bulk of a galaxy's flux. In addition, due to the convolution with a sub-pixel 
sampled PSF, and using the weighting scheme discussed above, we can reliably determine 
basic morphological parameters for even galaxies with strong quasar-like nuclei. 

5.3. Consistency 

We checked the consistency of the 1-D and 2-D fits by comparing the median ellipticity 
and position angle determined for the ellipse fits with the best-fit ellipticity and position angle 
in the parametric fits. For this comparison we exclude objects that have a significant nuclear 
point source, defined here as being one that contributes > 10% of the total flux of the source. 
The median position angles of each source can clearly be seen to be consistent with the fitted 
parametric values (Fig. H]), with three outliers: 3C 129, 3C 192 and 3C 321. The ellipticities 
show somewhat more variation between the two modeling techniques, with fourteen sources 
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differing by > 50%: 3C 79, 3C 88, 3C 105, 3C 130, 3C 135, 3C 192, 3C 264, 3C 284, 3C 319, 
3C 321, 3C 338, 3C 401, 3C 438, 3C 465 However, the Sersic parameters obtained from tfie 
2-dimensional parametric fit, and tfie 1-dimensional fit to the elhptical isophote intensities 
exhibit a correlation, but are seen to differ significantly in many individual cases. The vast 
majority of objects have Sersic indices measured in each way falling within a factor of 1.5 of 
each other. Three significant outliers have higher 2-D Sersic indices than 1-D: 3C 88, 3C 296 
and 3C 338. 

Ellipse fits allow a more detailed tracing of the structure, and thus have significantly 
better residuals than the 2-D Sersic model fits, which cannot cope with the twists in isophotes 
seen in many sources. Nevertheless, both techniques do a good job of modeling the bulk 
of the source flux in the majority of cases. We consider the 2-D radial profile fits to be 
the "correct" approach to obtaining a meaningful radial profile. The 1-D fit to elliptical 
isophotal data does not constitute a true radial profile due to overlapping adjacent elliptical 
annuli, and thus the form of the Sersic law determined from this technique can only be 
an approximation. Fitting with two or three Sersic components was attempted in several 
cases, and while this results in significant impro vement in the fi t , it is at the cost of physical 



interpretation of the resulting parameters. See iDonzelli et al.l (120071 ) for 2-component fits 
to the elliptical isophotes for a large number of objects in the present sample. Note that 
the introduction of a second fitted component can drastically alter the form (Sersic index, 
scale length, luminosity) of the primary modeled component. Donzelli et al. likely better 
model the true luminosity of the galaxy as their model provides a better overall fit to its 
form. However, they do not provide fits to the nuclei, and thus the quality of the Donzelli 
fit to the inner 1.7"is generally poor. The present technique provides a reasonable estimate 
of the luminosity and form of the bulge and of the brightest nuclei, and allows meaningful 
comparison with existing samples of quasar host galaxies and quiescent ellipticals in the 
literature. 

If we exclude all the nuclear-dominated sources (for which one cannot trust the ID fit 
results due to the asymmetric diffraction spikes dominating the ellipse-fitted isophotes) we 
find that the ID and 2D fitting results are essentially identical. Indeed, in simulations where 
we began with simulated galaxies and progressively added stronger and stronger nuclear 
components and random noise, we found that the 2D approach was always able to return the 
original parameters of the underlying galaxy, within the errors, even once the ID approach 
became unworkable due to the strength of the diffraction spikes, and resulting asymmetry in 
the isophotes. Furthermore, adopting only the non-nucleated sample which one can study 
using both techniques, we find no significant change in our overall results. 
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5.4. Host galaxy properties 

Fig. [5] sliows tlie Re — fie distribution resulting from each modeling technique, excluding 
two outliers, 3C 71 and 3C 258. 3C 71 has an unreliable fit due to its extension across both 
the N1C2 and WFPC2 fields (section S^H]). We beheve the highly pecuhar source 3C 258 to 
be a higher redshift quasar behind a. z = 0.165 irregular galaxy (see section 14.1.91) . Using the 
2D results for the reasons discussed above, the mean effective radius, or scale-length (defined 
throughout as the half light radius of the fitted models) for the entire sample is found to be 
(Re) = (-R1/2) = 7.46 ± 2.00 kpc (median= 3.63 kpc), with the mean surface within the half 
light radius, (/i)e = 18.12±0.25 mag.arcsec"^ (median= 17.59 mag.arcsec"^). Quoted errors 
are standard errors on the mean throughout. The mean host luminosity for the sample is 
{Mh) = —23.53 ± 0.27 mag (median=-23.75). The mean (2D) Sersic index for the sample 
is (n) = 3.68 ± 0.19 (median=3.41), and the mean ellipticity using the same technique is 
(1 - b/a) = 0.20 ± 0.02 (median=0.18). 

Excluding 3C 258 and 3C 71 we obtain the following 3-cr rejected means and standard 
errors for the sample: 

(Re) = 7.59 ± 2.00 kpc (median = 3.66 kpc). 

(/i)e = 17.89 ± 0.12 mag.arcsec"^ (median = 17.57). 

(Mh) = -23.81 ± 0.10 mag (median = -23.77). 

(n) = 3.69 ± 0.19 (median = 3.47). 

(1 - b/a) = 0.20 ± 0.02 (median = 0.18). 

We now briefly consider the main subdivisions of the sample into PR I and FR II types. 

5.4.1. FRI's 

For the FR I subsample, the 3-sigma rejected statistics cirG clS follows: 

{Re) = 12.03 ± 7.66 kpc (median = 2.74). 

(/i)e = 17.75 ± 0.30 mag.arcsec"^ (median = 17.49). 

(Mh) = -23.47 ± 0.27 mag (median = -23.29). 

{n) = 2.77 ±0.31 (median = 2.39). 

(1 - b/a) = 0.18 ± 0.02 (median = 0.18). 



-30- 




Mectiaa ellipUoity (ellipse fli) Median PA (elUpse Os) SerjiSe 



Fig. 4. — Results of the Ellipse (ID) and Galfit (2D) fits compared for targets without strong 
nuclear point sources (i.e. in which the nucleus contributes < 10% of the total luminosity). 
Left-to-right: EUipticity, position angle (PA), and Sersic index. We plot the median ellipticity 
and PA of the eUipse fits, which are found to be consistent with the 2D parameters. The 
Sersic profile indices, n, fitted to the elliptical isophotal intensities differ from those fitted 
directly to the data through the 2-dimensional parametric model by up to a factor of 1.5 
(indicated by the dotted lines). The 3 outhers are 3C 88, 3C 296 and 3C 338. 




Fig. 5. — Scale-length - Surface-brightness relation for the ID (left) and 2D (right) model 
fits to the full sample, excluding the outhers 3C 71 and 3C 258 (not plotted - see text). 
Symbols are the same as in Fig. 1. 
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The FR I's as a class tend to have lower n ("diskier") host galaxies than the combined 
sample, or any other individual subsample. They account for most of the lowest redshift 
[z < 0.05) sources, with few at higher redshifts. They exhibit a large spread in scale length, 
being found both in normal elliptical galaxies, and in cD galaxies like 3C 84. 

5.4.2. FR IPs 

For the FR II subsample the 3-sigma rejected statistics are as follows: 

{Re) = 5.94 ± 0.80 kpc (median = 3.84). 

(/i)e = 17.95 ± 0.13 mag.arcsec"^ (median = 17.61) 

(Mh) = -23.88 ± 0.10 mag (median= -23.84). 

(n) = 3.88 ± 0.21 (median = 3.70) 

(1 - b/a) = 0.21 ± 0.02 (median = 0.19). 

The FR II's represent a slightly more homogeneous population of giant "true" elliptical 
galaxies, with Sersic indices close to 4. The median size is larger than that of the FR I's, and 
the scale-length spread is much smaller. However, a two-sided Kolmogorov-Smirnov (K-S) 
test reveals little significant difference between the distributions in host luminosity or scale 
length. There is an 18% chance of the FR I and FR II Sersic indices being drawn from the 
same distribution {D = 0.26). 



6. Companions, jets, mergers and interactions 



We used the Source Extractor software ( Bertin and Arnoutslll996 ) to detect and measure 
each candidate companion source from the ellipse galaxy modeling residual images. We also 
obtained WFPC2 optical images of all of our sources from Andre Mart el's 3CR databasjil, o r 
directly frorn theMASlj^, most of which were previously published by de Koff et al. ( 1996 ): 



Martel et al.l (119991 ). Table [6] summarizes the origins of the optical images used in this 
study. We present basic companion astrometry and iJ-band photometry in table U\ along 
with an inventory of the number of companions to each 3CR host galaxy visible on the 



^http: / /acs. pha.jhu.edu/~niartel/ 
^http: / / archive.stsci.edu/ 
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NIC2 chip. We also manually investigated each object in the optical an d the infrared, to 
classi fy each according to the presence of the following features: Shells (IMalin and Carter 



19831 ): tidal tails; mergers (clearly undergoing merger with a companion from distortion or 



overlapping isophotes); pre-mergers (appear to be close to merging with a companion source); 
major companion sources (with infrared luminosity < 1 H mag. dimmer than the primary 
host galaxy assuming the same redshift); minor companion sources (a resolved companion 
source > 1 H mag. dimmer than the primary host galaxy assuming the same redshift); 
unresolved point-like companion source (a candidate companion source that is unresolved 
in the NICMOS image); jets, and dust disks. The table lists the presence of each type of 
artifact. Below we discuss the main types encountered. 



6.1. Globular Clusters 

In the very nearest sources, the large number of detected faint, unresolved sources are 
likely to be Globular Clusters within the halo of the host galaxy. This is really only an issue 
in the most well-studied, nearby objects, i.e. those from archival data. In particular, 3C 71 
(better known as NGC 1068) is well known for its globular cluster population. Other objects 
with candidate globular cluster detections are: 3C 29 {z = 0.05); 3C 31 {z = 0.02); 3C 84 
(Per A / NGC 1275; z = 0.02); 3C 88 {z = 0.03); 3C 129 and 3C 129.1 (both at z = 0.02). 
At present the sample is too small, but it would be interesting from the point-of-view of 
AGN feedback models to explore whether radio jets affect the formation of globular clus- 
ters, by examining the color distributions of the globular cluster populations in these objects 
compared with those in quiescent ellipticals. In the merger scenario for the formation of el- 
liptical galaxies, an older stellar population combining with a younger one produces a spread 
in metall icity and a birnodal d istribution in the metallicity and hence color of the globular 



clusters ( Ivan den Berghl |2000| ) . It is important and interesting to determine whether the 
globular cluster populations of radio galaxies differ from those of normal quiescent elliptical 
galaxies, but this lies outside the scope of the present study. 



6.2. Jets and Hotspots 



Two new opt ical-IR jets have be en discovered by the present NICMO S imaging survey; 
Those of 3C 133 Jpiovd et al.ll2006ah and 3C 401 Jchiaberge et al.ll2005h . 3C 130 exhibits 
numerous point sources, but is too distant (at z = 0.109) for them to be explained by 
Globular Clusters. Many of these may be hotspots, corres ponding to the u nusual radio 
morphology, and numerous hotspots observed in the radio by (lHardcastlelll998l ) . We identify 
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a further six sources with candidate IR-synchrotron hotspots, based on their distance, the 
number of unresolved sources, unusual radio morphologies and/or spatial coincidence of the 
unresolved sources with features on the radio map: 3C 52; 3C 61.1 (single bright unresolved 
source to north of galaxy); 3C 66B; 3C 219; 3C219.1 and 3C 452. 



^Jets are iden tifiable in a further ten sources: 3C 15 (jMartel et al.lll998l) : 3C 66B (IButcher et al, 

1980h : 3C 133 Jpiovd et all [2006a[): 3 C 264 JCrarie et all Il99'3f)~ 3 C 274 JCurtisI Il917h: 
3C 2 77.3 jMilev et al.lll98lh: 3C 293 Jpiovd et al.ll2006bh: 3C 346 JPev and van Breueel 
iggj; 3C 371 JmIssou et al.lll997h : and 3C 401 Jchiaberee e7aDl2005h Extended emission 



line regions, with the appearance of jets are clearly visible in 3C 171 and 3C 234. 



6.3. Dust disks 



31 objects in the sample are seen to possess ~ 100 pc scale distributions of dust and 
gas in their nuclea r regio ns, often in settled, disk-like structures surrounding the AGN. 
See iTremblay et al.l (120071 ) for further details. 11 of these are FR I sources; 16 are identified 
as FR IPs. 



6.4. Mergers and pre-merger candidates 

Six sources are clearly identifiable as mergers through their overlapping isophotes and 



disturbed appearance: 3C 79; 3C 293 flFlovd et al.l l2006af ): 3C 321; 3C 346; 3C 405 and 



3C 433. There are an additional twenty-eight sources that qualify as "pre-mergers", or can- 
didate mergers that cannot be confirmed by the overlapping of isophotes of the two partners. 
"Major" secondary sources are identifiable in twenty sources (20%). Minor secondaries are 
identified in fifty-five sources, or 56% of the sample. Fifty-nine sources have unresolved 
companions. A number of these may turn out to be foreground stars, especially for sources 
at low galactic latitudes, but it seems unlikely that all such sources are so explained. For 
a small number we also obtain detections in the optical, and these sources typically have 
colors of R — H = 2 — 3. The exact nature, and redshifts of these sources must await a 
detailed spectroscopic study, which could confirm whether the candidate companions indeed 
lie at the same redshift as the radio source, and further identify the stellar makeup of the 
sources. However, based on the sheer number of small companions, we conclude that these 
are of significant importance to the 3CR sample and to the radio source phenomenon itself. 
The broad-band colors of many of these sources are consistent with an old stellar popula- 
tion, perhaps the stripped core of a late-type spiral, or compact early-type galaxy that has 
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undergone merger with the main galaxy, losing its gas and dust to the more massive system. 
Such a mechanism would provide a natural means to fuel the AGN activity. As far as we are 
aware, this is the first such detection in large numbers of sources across a range in redshifts. 
Examples of dwarf ellipticals and of globular clusters have been observed in a number of 
nearby sources, e.g. Fornax. 

7. Discussion 

7.1. Galaxy morphologies 

We find a broad distribution of Sersic indices in our sample, spanning the full range 
from n = 1 (exponential disks) to n = 4 (de Vaucouleurs ellipticals), and higher (see Fig. [H]). 
There is a moderately strong correlation (Spearman's rank p = 0.42, p = 1.7E — 5) between 
scale- length and Sersic n (see Fig. [7]). Interestingly, neither the FR I's nor the FR II's 
considered alone show such a significant correlation (FR I: p = 0.56, p = 0.0055; FR II: 
p = 0.35, p = O.OOTS-E — 3). There is a somewhat weaker anti-correlation between the host 
galaxy absolute magnitude and the Sersic index, n (p = — 0.32,j9 = 0.001) - once again 
somewhat weaker in the individual subsamples. Thus the larger, more luminous galaxies 
tend to be more "bulgy". The Sersic index is also found to correlate weakly with redshift 
(p = 0.28; p = 0.006). Thus the disky sources tend to be lower redshift, and lower IR 
luminosity, but interestingly, no strong relationship is found between Sersic index and radio 
power (either core or total). However, there are very few sources that have Sersic indices 
n < 2, thus this is not an indication that true "disk" galaxies are capable of hosting powerful 
radio sources. Rather, it indicates that within the radio source population, there is a strong 
requirement for an elliptical, or bulge-dominated host galaxy, and a strong dependence of 
radio luminosity on host galaxy size and luminosity, but not on Sersic index within the range 
2 < n < 10. A more "bulgy" bulge does not allow for a more powerful radio source, while a 
more luminous bulge does. 

There are a number of exceptions to the general rule; Seventeen of the host galaxies have 
n < 2, and would thus be classed as disky. Ten of these sources are low-redshift {z < 0.05), 
nine of those being FR I low-power radio sources. However, there are eight FR II sources, 
and five of these are at ^ > 0.1: 3C 173.1, 3C 234, 3C 288, 3C 323.1 and 3C 332. All five of 
these FR II sources have close companion objects, and are either merging, or are likely to be 
post-merger systems. They generally show evidence of dusty environments from their optical 
and IR images, except for 3C 323.1 which has a quasar-like nucleus in both the optical and 
IR. 
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Fig. 6. — Host galaxy morphology distribution in terms of Sersic index, n: FR I sources are 
shown shaded, with the sohd hne indicating the entire sample. The host galaxies exhibit a 
broader range of Sersic n than a population of pure ellipticals, with a tail of low-n galaxies 
(nineteen at n < 2). These "disky" sources are hsted and briefly discussed in the main text. 



- 36 - 



7.1.1. Ellipticity 

The 3CR host galaxies are shghtly rounder, on average, than the general elliptical galaxy 
population, with a sharp peak at El, and no objects more eccentric than E5. See Fig. [HI 



7.1.2. Isophotal twists 

The total isophotal twist angle, was computed for each galaxy from the position 
angle profile. We neglect the inner Of! 5, and neglect > la outliers in the 6 profile. We then 
quantify the twist in terms of the total position angle variation of the major axis over the 
entire range in surface brightness (after iGopal-Krishna et al.ll2003l ). Objects that are less 
elliptical than 1 — b/a < 0.2 over their entire range are also discounted, as for these extremely 
round objects, a spuriously large twist angle can arise due to the circular symmetry of the 
isophotes. The resulting twist angles are plotted in Fig. [9l We confirm that th e trend 
between total radio po wer and twist angle f irst noted by iGopal-Krishna et al.l (120031 ) for the 
radio galaxy sample of iGovoni et al.l (120001 ) is also seen in the present sample. 



7.2. Galaxy luminosities and scale lengths 

The host galaxy luminosity distribution is shown in Fig. [TD], with their scale length 
distribution shown in Fig. [TTl All of the 3CR sources studied, except for 3C 258, 3C 270, 
and 3C 296, are hosted by galaxies with luminosities above Mh = —22, and scale lengths 
> 1 kpc. The highly peculiar and disturbed 3C 258 (see earlier comment in section 14.1.91) 
is found to have an extremely compact and low-luminosity host galaxy, consistent with it 
being at a much higher redshift. 

Half of the sources have host galaxies brighter than L*, and larger than 3 kpc, with the 
peak in the luminosity and scale-length distributions just below L* an d 3 kpc respectively . 



This is as to be expe cted based on our knowledge of radio galaxies (e.g. IZibetti et al 



20021 



i 



quasa r host galaxies ( Dunlop et al. 2003 : Floyd et al. 2004) and unification schemes (jSarthel 
19891 ): Our understanding of the AGN phenomenon hinges upon the presence of a supermas- 
sive black-hole as the ultimate energy source, and our growing demographic understanding 
of galaxies is that a massive black-hole requires a massive galaxy bulge or spheroid to host 
it. 



The low-luminosity tail of the sample is dominated by low-redshift and low-radio-power 
FR I's: Of the twenty-two sources with luminosities Mh > —23, nineteen are at z < 0.09, 
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with only 3C 61.1, 3C 258 (once again), and 3C 314.1 breaking the general rule. 



3C 61.1 at z = 0.186 has an unusual radio morphology, and appears to be hosted by 
an unusually faint, small galaxy at th e center of a small group. The optical image shows 
complex structure (jde Koff et al.lll996l ). with tails of emi ssion that resemble spira l arms, and 
which do not correlate with the radio emission features (ILeahy and Perleylll99ll ). However, 
the published 20 cm VLA radio map is of low resolution, and it is difficult to confirm the 
location of the optical-IR counterpart accurately. 



3C 314.1 has an unusually low-surface-brightness host galaxy at z 
slightly dusty appearance from the optical-IR images. 



0.119, with a 



7.3. Infrared nuclei 



A number of our sources (40%) have clearly detectable infrared nuclear activity. We 
were able to detect nuclear point sources down to an absolute magnitude of Mh{Nuc) = — 15 
(table O Fig. [T2|) . A more exhaustive search for faint nuclei, and a discussion of their origin 
will be the subject of another paper (Chiaberge et al. in preparation). Firstly, we note that 
nuclear point sources are detected in both FR I's and FR IPs, with roughly equal probability 
within our sample (38% of FR I's, 41% of FR IPs). No significant correlation is found 
between the NIR nuclear luminosities of the sources in which we have identified unresolved 
nuclear point sources, and the host luminosity, nor scale-length, although a weak correlation 
is seen with host galaxy surface brightness (p = 0.34; p = 0.035) - see Fig. [131 The difference 
in the luminosity range of FR I and FR II sources is very clearly seen from Fig. [T21 FR I's 
although they can possess nuclei, do not contain luminous quasar-like nuclei that are seen 
in the FR II population. 



7.4. Kormendy relation (R-/i) 



The strong correlation between scale len gth and surface brightness is a well-known 
feature of elliptical galaxies (IKormendyl 119771 ) . The relation followed by the present low- 
redshift 3CR sample is somewhat lower tha n that of 3.33 ± 0.09 established for early type 
galaxies using the Sloan Digital Sky Survey (IBernardi et al.ll2003l ). following a slope of 3.0 ± 
0.1. 



/i3CR = 16.0 ± 0.1 + (3.0 ± 0.1) logio(i?e) 
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Kormendy fit to the FR I's and FR II's yield somewhat different slopes of 3.46 ± 0.1 
and 2.67 ±0.1 respectively (Fig. [T^ . 

fiFRi = 15.9 ± 0.1 + (3.46 ± 0.1) logio(i?e) 

(IFRH = 16.1 ± 0.1 + (2.67 ± 0.1) logio(i?e) 



7.5. Fundamental plane 

We used Principal Component Analysis (PCA) to explore the model parameter space 
for any 3-dimensional correl ations, like the photomet ric "fundamental plane" reported for 
normal elliptical galaxies by ( Khosroshahi et al. 2000bl lal). PCA looks for correlations among 



a set of multi-dimensional data and rotates the coordinate axes from being aligned with the 
input parameters to new directions that are aligned with the directions of greatest variance of 
the data-the principal axes. In order to prevent the choice of units from artificially weighting 
some parameters more than others, each parameter is first normalized by subtracting its mean 
and dividing by its variance. The PCA is then performed on these normalized variables. In 
the output, each of the eigenvectors - the principal axes- is written as a linear combination 
of the original (but normalized) parameters. The eigenvalues are scaled so that the sum of 
all eigenvalues equals the total number of eigenvectors (and therefore the total number of 
parameters as well). In Fig. [TS] we show the "scree plot" of eigenvalues associated with the 
seven eigenvectors fit, along with a plot of the data over the first eigenvector. Only the first 
eigenvector accounts for a significant quantity of the scatter in the parameter space, and is 
made up almost entirely of Re, f^e, and Sersic n. The recovered form of the relationship is: 

= 2.771ogio Re + 1.731ogio n + 14.78 



7.6. Host galaxy colors 

For the cross-sample of this paper and Martel et al. (1999) we can produce accurate 
R — H colors for the host galaxy, shown in table [81 A color - absolute magnitude diagram 
for this subsample is presented in Fig. [161 An extremely steep trend is seen between color 
of the host galaxy {R — H) and absolute if-band magnitude of the host galaxy, with a slope 
of roughly 1. The extremely red colors at the top end of the sequence are likely due to dust 
absorption in the optical images of sources like 3C 293, 3C 296, 3C 388, and 3C 403. The 
blue outliers are all nucleated sources, in which the Martel et al. fluxes are likely to be higher 
than the true galaxy flux, as no attempt was made by those authors to separate host and 
nucleus. 
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7.7. Host galaxy masses 

Based on the apparent similarity between the present sample, and the nor mal quiescent 



elliptical galaxy population, we have applied the finding of IZibetti et al.l (120021 ) that if-band 
luminosity provides a first order estimate of the dynamical mass of an elliptical galaxy. They 
find a mean mass-to-light ratio of log(M/L/^) = 0.09 (in solar units), that is independent 
of total luminosity. There is an approximately 0.2 dex la scatter on this relationship. The 
mean mass of a 3CR host galaxy from the present sample is 2 x 10^^ Mq, but with very 
large scatter, in particular due to the 2 enormous outliers, 3C 130 (10^^ ^o) and 3C 338 
(5 X 10^^ Mq). Clearly the assumption of constancy of the mass-to-light ratio for our sample 
needs verification by dynamical means, and we would argue that a significant hole in the 
current 3C literature is a fundamental plane relationship (R-fi-cr) study to test whether these 
objects truly are dynamically representative of normal elliptical galaxies. In the absence of 
long- slit spectra in the public domain such a work is currently non-existent. 



8. Summary and conclusions 

Ellipse and Galfit models have been presented for 101 z < 0.3 3CR radio source host 
galaxies. The two fitting methods (ID and 2D) are found to be statistically consistent for 
non-nucleated sources, although individual sources can provide significantly differing results 
using the two methods. Simulations have revealed that, with the addition of an artificial 
nuclear point source to a real source, we are still able to obtain accurate results using the 
2-D technique, after the 1-D technique has failed. 

Unresolved nuclear point sources are detected with equal likelihood down to Mh = —15 
in both FR I and FR II galaxies. FR I's can have faint unresolved nuclei (all dimmer than 
Mh = —20), but these mostly exist in low-redshift sources in this radio-fiux limited sample. 
We find no significant difference in properties of the host galaxies between nucleated and 
non-nucleated sources. 

The host galaxies of the 3CR are generally consistent with the elliptical galaxy pop- 
ulation: They exhibit a Kormendy relation that is similar to that of quiescent ellipticals 
and quasar host galaxies; the peak in the host galaxy luminosity distribution is close to 
L*, with steeper drop-off to high luminosity than low; there is a strong correlation between 
host if-band luminosity and {R — H) color. However, the 3CR exhibit a large spread in 
Sersic index to low values, and includes a number of quite disky sources. Sersic n is not 
found to correlate with the radio luminosity, and yet the disky (low-n) wing of the sample is 
dominated by low-2, low radio-power sources, merging sources, and sources with luminous 
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companions. This suggests that disky galaxies are capable of hosting powerful radio sources, 
providing that they are massive enough, or are undergoing a major merger. 

Approximately 50% of our sample exhibit compact, often unresolved, companion sources 
that are too red and faint to have been detected in previous optical snapshot programs. 
Spectroscopy is required to determine the nature and redshifts of these sources. Their typical 
colors (where this can be determined) are around R—H = 2—3, consistent with mature stellar 
populations at the redshift of the primary source. These sources are found predominantly 
in elliptical host galaxies, and have 3 likely origins: Foreground stars - likely in cases with 
low galactic latitude; Synchrotron Hotspots - particularly in a few peculiar radio sources 



Merger remnants - compact cores from cannibalized s mall galaxies (e.g. ICanalizo et al. 



20031 ) : Or molecular gas clouds infalling into the galaxy (iBellamy and Tadhunterll2004l ). By 
comparison, only ~10% of our sample show signs of an ongoing or recent major merger, 
typical of the elliptical galaxy population in general (though many more show other signs of 
disturbance). If many of these sources do turn out to be galactic nuclei, it would suggest that 
a minor merger is sufficient to fuel or re-fuel a quiescent black hole in an elliptical galaxy into 
radio-loud AGN activity, while disky galaxies appear to require a more major disturbance. 
There is a correlation between radio power, and the isophotal twist angle within our sample, 
with all sources exhibiting significant isophotal twists. We thus conclude that while bulgy 
elliptical galaxies host the majority of powerful radio sources, it is possible to trigger such 
activity in diskier objects through mergers. The diskiness may simply be a morphological 
distortion resulting from the merger, or the merger may be able to provide sufficient fuel to 
a black-hole in a disky galaxy to trigger behavior only normally seen in the more massive 
black holes at the centers of ellipticals. 

A detailed spectroscopic study of the 3CR host galaxies is overdue. It is important 
to explore the dynamical states of these objects to place them on the "fundamental plane" 
{R — ji — a) of galaxies. The nature of the companion sources needs to be explored, ideally 
using imaging spectroscopy to determine their redshifts and compositions at the same time 
as exploring the makeup of the 3CR host galaxy itself. 
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Fig. 7. — Weak correlations are observed between the Sersic index, n, and the host galaxies' 
luminosity and scale length. Powerful radio galaxies can only be hosted by disk- dominated 
galaxies if the galaxy is exceptionally massive. Symbols are the same as in Fig. 1. 



Table 1. Observing log 



Source 






■^2000 




z 


X 


Obs. Uatc 


3C 


15 


00 


37 


04 


1 


-01:09:08 





073 


-64 


2006-06-24 


3C 


17 


00 


38 


20 


5 


-02:07:40 





219 


-65 


2006-07-02 


3C 


33 


01 


08 


53 


3 


-1-13:20:25 





059 


-49 


2006-06-29 


3C 


98 


03 


58 


54 


4 


-H10:26:03 





030 


-31 


2005-11-25 


3C 


132 


04 


66 


43 





-H22:49:22 





214 


-13 


2005-11-24 


3C 


153 


06 


09 


32 


6 


-H48:04:16 





277 


13 


2005-11-25 


3C 


166 


06 


46 


24 


1 


-1-21:21:61 





245 


8 


2005-11-04 


3C 


234 


10 


01 


49 


5 


-1-28:47:09 





185 


53 


2005-11-03 


3C 


258 


11 


24 


43 


5 


4-19:19:12 





165 


69 


2006-01-25 


3C 


284 


13 


11 


04 


7 


-1-27:28:08 





239 


86 


2006-03-04 


3C 


296 


14 


16 


52 


9 


-1-10:48:26 





025 


64 


2006-04-21 


3C 


300 


14 


23 


01 





-M9:35:17 





270 


68 


2006-03-04 


3C 


323.1 


15 


47 


43 




4-20:52:17 





264 


49 


2006-04-21 


3C 


326 


15 


52 


09 


i 


-1-20:05:24 





090 


48 


2006-04-21 


3C 


332 


14 


23 


01 





4-19:35:17 





270 


45 


2006-01-12 


3C 


357 


17 


28 


18 


5 


4-31:46:14 





166 


31 


2006-03-25 


3C 


403.1 


19 


52 


30 


4 


-01:17:20 





055 


-14 


2006-06-23 


3C 


410 


20 


20 


06 


6 


4-29:42:14 





248 


-4 


2006-06-30 


3C 


424 


20 


48 


12 


1 


4-07:01:18 





127 


-22 


2006-06-22 


3C 


442 


22 


14 


46 


9 


4-13:50:27 





026 


-34 


2006-06-27 


3C 


459 


23 


16 


36 


1 


4-04:06:19 





219 


-51 


2006-06-24 



Note. — Observation dates and positions of the 21 newly observed 
sources in the sample. We present the source name, J2000 Equatorial 
coordinates, redshift, galactic latitude, and observing date. Two other 
targets, 3C 18 and 3C 63 were observed as part of the program but 
were missed by the NIC2 chip. 
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Fig. 8. — EUipticity distribution for the full sample (unshaded) and for the FR Fs (shaded). 
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Fig. 9. — Isophotal twist angle, A6' against (left-to-right): total radio power, boxiness (04), 
and median ellipticity. Round er galaxies tend to have larger isophotal twists, and as found 
by (IGopal-Krishna et al.ll2003l ). radio galaxies exhibit a correlation between twist angle and 
total radio power. Symbols are the same as in Fig. 1. 



Table 2. Archival Objects: Observation Dates 



Source 


Alt. names 


"2000 


•^2000 


2 


A 


Obs. 


Datc 


PROPOSID 


3C 71 


M 77, NGC 1068 


02:42:40.7 


-00:00:48 


0.003793 


-52 


1998- 


02-21 


7215<"' 


3C 84 


Per A, NGC 1275, A 0426 


03:19:48.1 


+41:30:42 


0.017559 


-13 


1998- 


03-16 


7330<''' 


3C 264 


NGC 3862, A 1367 


11:45:05.0 


+ 19:36:23 


0.021718 


73 


1998- 


05-12 


7862 


3C 270 


NGC 4261 


12:19:23.2 


+ 05:49:31 


0.007465 


67 


1998- 


04-23 


7868''** 


3C 272.1 


M 84, NGC 4374 


12:25:03.7 


+ 12:53:13 


0.003536 


74 


1998- 


07-13 


7868'=' 


3C 274 


M 87, NGC 4486 


12:30:49.4 


+ 12:23:28 


0.004360 


75 


1997- 


11-20 


717l(/) 


3C 293 


UGC 08782 


13:52:17.8 


+31:26:46 


0.045034 


76 


1998- 


08-19 


7853(9' 


3C 305 


UGC 09553, IC 1065 


14:49:21.6 


+ 63:16:14 


0.041639 


49 


1998- 


07-19 


7853'''' 


3C 317 


UGC 09799, A 2052 


15:16:44.5 


+ 07:01:17 


0.034457 


50 


1998- 


08-26 


7886'''' 


3C 338 


NGC 6166, A 2199 


16:28:38.5 


+39:33:06 


0.030354 


44 


1997- 


12-18 


7453('' 


3C 405 


Cyg A 


19:59:28.3 


+40:44:02 


0.056075 


6 


1997- 


12-16 


7258'^' 



Note. — Observation dates and positions of the 11 archival objects in the sample. All observations are on NIC2 
through the F160W filter. We list HST proposal ID's and references for the original publication of the data. 3C 273 
also has deep NICMOS F160W imaging of its jet, but is excluded from this paper as the on-galaxy integration time is 
too short for an accurate characterisation of the host galaxy. 



References. — (a) [Thompson and CorbinI l |1999|): (b) iMartini et al.| <2003|): ( c ) iCapetti et al.| l|2000bD: 

fd) | QuiIIen et al.| l|2000r }; (e) [ Bower et al.| i200GD: ff) JCorbin et al.| 1 I2OO2I '}; fg) I Floyd et al.| l |2006bD : fh) I Jackson et al.| 
12003) : fi) lJensen et^ i200l|) : fi ) iTadhunter et al.| l |1999l '). 
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Fig. 10. — Galaxy luminosity distribution for the full sample (unshaded) and for the FR I's 
(shaded): SCR sources generally hosted by super-L* objects. The low- luminosity tail is 
dominated by low-redshift {z < 0.09) and FR I targets. 
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Fig. 11. — Scale-length distribution for the full sample (unshaded) and for the FR Fs 
(shaded) . 
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Fig. 12. — Nuclear luminosity distribution for the full sample (unshaded) and for the FR Fs 
(shaded): We have detected unresolved IR nuclei in 39 objects, down to an absolute H-band 
magnitude of —15. 
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Fig. 13. — Nuclear IR luminosity vs. host galaxy properties, left-to-right: Scale-length, 
surface brightness, and host galaxy luminosity. There is a slight tendency for stronger 
nuclear sources to be located in higher surface-brightness galaxies. No trend is observed 
between nuclear luminosity and host luminosity (right) or scale-length (left). 



Table 3. Radio Properties of the NICMOS Snapshot Survey Sample 



Source 


2 


V 


S(178 MHz)"/Jy 


logio Pits 




las''-'' 






Olass'' 


^OII 


^01 II 


3C 15 


0.073 


15.3 


15.8 


26.3 


0.64 














3C 17 


0.219 


18.0 


20.0 


27.4 


0.52 














3C 20 


0.174 


19.0 


42.9 


27.6 


0.67 


51 


101 


"2 


HEG 


33.99 


34.47 


30 28 


0.195 


17.6 


16.3 


27.2 


1.06 


30 


166 


2 


LEG 


35.07 




3C 29 


0.044 


14.1 


15.1 


25.8 


0.50 


139 


160 


1 








3C 31 


0.016 


12.1 


16.8 


25.0 


0.57 


1833 


159 


1 


FRI 






3C 33 


0.059 


15.2 


54.4 


26.7 


0.76 






2 


HEG 


34.72 


35.37 


3C 33.1 


0.181 


19.5 


13.0 


27.1 


0.62 


2i6 


45 


2 


WQ 




35.05 


3C 35 


0.067 


15.6 


10.5 


26.1 


0.77 


704 


12 


2 


LEG 






3C 52 


0.285 


18.5 


13.5 


27.5 


0.62 


51 


20 


2 








3C 61.1 


0.186 


19.0 


31.2 


27.5 


0.77 


186 


2 


2 


HEG 


34.70 


35.48 


30 66B 


0.021 


12.9 


24.6 


25.4 


0.62 


330 


54 


1 


FRI 


33.16 


33.01 


30 71 


0.004 


8.9 


16.1 


23.7 


0.55 














30 75N 


0.024 


13.6 


25.8 


25.5 


0.71 


692 


lii 


i 








30 76.1 


0.032 


14.9 


12.2 


25.5 


0.77 


200 




1 


FRI 







Note. — Main radio and optical properties for each source: rcdshift; V-band magnitude; radio flux density in Janskys at 
178 MHz; radio power at 178 MHz in W Hz~"''; radio spectral index; radio source largest angular size; radio structure position 
angle; Faranoff-Riley class; Ionization Class; [OH] 3727 emission line luminosity in W; [OIII] 5007 emission line luminosity in W. 
[The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample.] 



Ref erences. — " ISpinrad et a l.|jj^9 85|) - updated v alues were taken from the NASA Extragalactic Database (NED); ^ |de Koff et all 
Jl996|) and references therein; '^" |A^i!rtel et all i 19991 ) and references therein; [Jackson and Rawlingj il997|) . 
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Fig. 14. — Best fitting scale- length-surface brightness relation for the full sample (solid line), 
the FR I's (dashed) and the FR IPs (dot-dash). The locus of an L* galaxy is shown by the 
dotted line. 
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Fig. 15. — Left: Eigenvalues from the Principal Component Analysis on the fitted parameters 
from the 2-dimensional (Galfit) modeling. Only the first eigenvector accounts for a significant 
proportion of the scatter in the entire dataset. Right: Photometric fundamental plane 
(eigenvector 1). 







Table 4. 


Ellipse and 1-D Sersic model fits 








Source 


(1 - b/a). 










Tiin '^4/*lma:r 


re /kpc 




We 


n 


3C 15 


0.06 


21.32 


-0.08 


0.02 


-0.04 


0.02 


4.96 ± 0.45 


19 


76 ± 0.13 


4.17 ± 0.08 


3C 17 


0.17 


21.32 


-0.17 


0.06 


-0.08 


0.14 


5.84 ± 1.53 


20 


69 ± 0.62 


10.00 ± 0.00 


3C 20 


0.10 


-63.77 


-0.06 


0.10 


-0.03 


0.20 


2.72 ± 1.45 


19 


17 ± 1.06 


3.23 ± 0.64 


3C 28 


0.18 


-26.43 


-0.03 


0.11 


-0.04 


0.03 


6.74 ± 0.63 


20 


26 ± 0.11 


2.01 ± 0.08 


3C 29 


0.03 


-76.28 


-0.05 


0.02 


-0.05 


0.03 


4.44 ± 0.15 


19 


05 ± 0.02 


1.50 ± 0.02 


3C 31 


0.12 


-32.23 


-0.02 


0.08 


-0.03 


0.01 


1.46 ± 0.05 


17 


20 ± 0.02 


1.43 ± 0.02 


3C 33 


0.07 


-43.84 


-0.09 


0.01 


-0.04 


0.06 


2.75 ± 0.21 


18 


97 ± 0.11 


3.19 ± 0.07 


3C 33.1 


0.06 


13.96 


-0.02 


0.03 


-0.05 


0.02 


2.19 ± 0.65 


18 


66 ± 0.65 


4.49 ± 0.33 


3C 35 


0.23 


-65.28 


-0.18 


0.02 


-0.07 


0.03 


4.59 ± 0.67 


19 


89 ± 0.20 


3.60 ± 0.14 


3C 52 


0.31 


10.96 


-0.03 


0.06 


-0.02 


0.08 


3.55 ± 0.36 


18 


74 ± 0.18 


1.73 ± 0.09 


3C 61.1 


0.09 


0.95 


-0.17 


0.22 


-0.07 


0.08 


1.87 ± 0.16 


19 


65 ±0.20 


2.07 ± 0.11 


3C 66B 


0.06 


-50.32 


-0.08 


0.04 


-0.03 


0.06 


2.42 ± 0.15 


18 


42 ± 0.04 


1.97 ± 0.05 


3C 71 


0.18 


36.54 


-0.04 


0.03 


-0.06 


0.02 


2.21 ± 0.51 


19 


51 ± 0.35 


10.00 ± 0.00 


3C 75N 


0.19 


-80.47 


-0.01 


0.05 


-0.02 


0.04 


1.64 ± 0.22 


17 


75 ± 0.16 


3.09 ± 0.12 


3C 76.1 


0.14 


-47.79 


-0.02 


0.04 


-0.05 


0.02 


1.12 ± 0.07 


17 


95 ± 0.09 


2.26 ± 0.06 


Note. - 


— Properties 


of the elliptical 


isophotal fit to each 


source, including 1-D radial-prohle Sersic fit; 


median ellipticity; median 



position angle (East of North); minimum value of 03 /a; maximum value of 03/0; minimum value of 04/0; maximum value of a^ja. scale 
length (half-light radius) of best-fit model to elliptical isophotes; mean surface brightness of model within scale-length, (/J.)e! Sersic index 
of model. [The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample.] 
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Fig. 16. — Color - absolute magnitude diagram for the sample. More luminous hosts are 
significantly redder. Nucleated sources (with nuclei brighter than Mh — —21) are shown by 
a filled star, and dominate the blue end off the distribution. 
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A. Consistency of Galfit and 2DM 



In this paper we liave used Galfit (jPeng et al.ll2002l ) to perform 2-dimensio nal modeling 
of the host galaxies of th e SCR. But we have adopted the error treatment of 2DM (IMcLure et al 
19991 : iFloyd et al.l |2004| ) in cases where the nucleus dominates (i.e. provides > 1/2 of the 
total flux). In this section we explore the consistency of the two approaches in a series of 
test cases using both real and simulated galaxies. 

2DM was designed for use on QSO host galaxies, and is a simpler fltting code than 
Galflt. The major difference is in assignment of errors, as for the PSF-dominated cores of 
QSO images it is the PSF sampling error that dominates any fltted model, and the Poissonian 
err or on the flux in a given pixel is an insuflicient determination of the error in the model 
flt (IFloyd et al.ll2004j ). In its standard form it flts a single Sersic host galaxy component and 
a nuclear component, centered at the same flxed position. 

Galflt and 2DM were used to model the host galaxy flux of three radio galaxies from the 
sample: 3C 66B, 3C 111, and 3C 449. These three were speciflcally chosen to span the range 
of infrared nuclear activity seen across the sample: 3C 449 shows no sign of an unresolved 
nuclear point source, while 3C 66B shows a weak point-like nucleus, and 3C 111 is a quasar. 
All data has been reduced following the technique described in section [21 We have modelled 
each object using a single Sersic component model in Galflt, as described in section [3], and 
using 2DM as described in lFloyd et al.l (120041 ). The resulting models are presented in table [91 
Next, synthetic galaxies were constructed based on the parameters recovered from these three 
objects, using Galflt. Random noise was added to each synthesized test case at incrementing 
levels, and these noisy synthetic targets are then remodeled with Galflt and 2DM to explore 
the stability of flt to noise, and ability to recover the "true" properties of the host galaxy. 
The results of this test are explored in section IA.2I below. 



A.l. Models of the three real test cases 

The best-flt model parameters for the three objects, 3C 66B, 3C 111, and 3C 449, using 
each of the modeling techniques, are presented side-by-side for comparison in table [3 Radial 
proflles are shown in Fig. [T71 2DM and Galflt are clearly seen from the table to be in good 
agreement. The host fluxes generally agree well (to within ~ 2%) and the scale-lengths 
within ~ 10%. However, in the case of the quasar, 3C 111, there is a larger discrepancy, 
and greater dependence on the treatment of errors in the central region. There is also a 
signiflcant discrepancy in the value of the nuclear flux in all three sources. This is once 
again, very dependent on the error treatment. 
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Table 5. 2-D Sersic model fits 



Source 


Mh-(Nuc) 


Mi, (Host) 






MH 










n 








1 










6 












d 






3C 


15 




-24.33 ± 0.01 


8.41 ± 


05 


18 


86 ± 





01 


6 


10 


± 





02 





07 


± 





00 


19 


33 


± 





38 





11 


± 





01 


3C 


17 


-21.51 


-24.19 ± 0.01 


4.25 ± 


05 


17 


52 ± 





00 


7 


49 


± 





11 





23 


± 





00 


20 


94 


± 





31 


-0 


21 


± 





02 


3C 


20 




-23.61 ± 0.00 


2.19 ± 


01 


16 


66 ± 





00 


2 


70 


± 





01 





10 


± 





00 


-65 


53 


± 





62 


-0 


15 


± 





01 


3C 


28 


-18.71 


-24.70 ± 0.01 


11.12 ± 


07 


19 


10 ± 





01 


2 


32 


± 





02 





24 


± 





00 


-33 


50 


± 





23 


-0 


26 


± 





01 


3C 


29 


-16.19 


-23.82 ± 0.01 


4.24 ± 


01 


17 


89 ± 





01 


1 


45 


± 





00 





04 


± 





00 


-81 


55 


± 





00 





08 


± 





00 


3C 


31 




-23.29 ± 0.00 


1.72 ± 


00 


16 


46 ± 





00 


1 


67 


± 





00 





12 


± 





00 


-32 


98 


± 





08 


-0 


01 


± 





00 


3C 


33 




-23.43 ± 0.00 


2.98 ± 


01 


17 


51 ± 





00 


3 


60 


± 





01 





08 


± 





00 


-35 


60 


± 





26 


-0 


13 


± 





01 


3C 


33.1 


-20.86 


-23.84 ± 0.02 


2.68 ± 


02 


16 


86 ± 





02 


6 


41 


± 





00 





01 


± 





01 


54 


91 


± 


2 


56 





00 


± 





00 


3C 


35 




-24.07 ± 0.00 


8.38 ± 


04 


19 


11 ± 





00 


4 


79 


± 





01 





28 


± 





00 


-66 


34 


± 





07 


-0 


02 


± 





00 


3C 


52 




-25.06 ± 0.00 


5.22 ± 


02 


17 


10 ± 





00 


2 


49 


± 





01 





35 


± 





00 


10 


16 


± 





18 


-0 


17 


± 





02 


3C 


61.1 


-16.96 


-22.18 ± 0.11 


1.73 ± 1 


56 


17 


59 ± 





00 


2 


83 


± 





00 





13 


± 





01 


-2 


02 


± 


1 


22 





41 


± 





00 


3C 


66B 




-23.35 ± 0.01 


2.93 ± 


01 


17 


56 ± 





01 


2 


39 


± 





01 





07 


± 





00 


-59 


24 


± 





34 





00 


± 





00 


3C 


71 




0.86 ± 0.01 


1.17 ± 


01 


39 


77 ± 





00 


3 


12 


± 





02 





38 


± 





00 


44 


99 


± 





09 





00 


± 





00 


3C 


75N 




-23.35 ± 0.00 


1.68 ± 


00 


16 


34 ± 





00 


3 


68 


± 





01 





14 


± 





00 


-78 


88 


± 





15 





22 


± 





01 


3C 


76.1 




-22.28 ± 0.00 


1.31 ± 


00 


16 


88 ± 





00 


2 


95 


± 





00 





16 


± 





00 


-48 


12 


± 





12 


-0 


08 


± 





00 



Note. — Properties of the best-fit 2-D Scrsic model to each source: Absolute magnitudes of the nuelcus and host galaxy; sealc-lcngth 
(half-light radius) in kpc Re', surface brightness at the scale- length fj,; Sersic index n; cllipticity 1 — — ; position angle 6 (relative to North); 
diskiness, d. diskiness, d. Notes: The following objects did not admit a complete or normal fit, as discussed in the main text: (a) 3C 71 
is too extended to be properly constrained on either the NIC2 or WFPC2 chips; (b) 3C 258 — values are presented as if the target lay at 
z — 0.165, however it is likely that the source is actually at higher rcdshift. [The complete version of this table is in the electronic edition 
of the Journal. The printed edition contains only a sample.] 



Table 6. Optical data from the WFPC2 archive 



Source 



PROPOtilU 



"Tilt; 



cr 



cxptimc 



3C 15 
3C 17 
3C 20 
3C 28 
3C 29 
3C 31 
3C 33 
3C 33.1 
3C 35 
3C 52 
3C 61.1 
3C 66B 
3C 71 
3C 75 
3C 76.1 



6348,5476 
6967,5476 
5476 
5476 

6967,5476 

6673,5476,6673 

5156 

6967,5476 

6967,5476 

6967,5476 

6348,5476 

6673,5476,6673 

5754,5754,5754,5479 

5476,5927 

6967 



F555W,F702W 
F555W,F702W 
F702W 
F702W 

F555W,F702W 

F555W,F702W,F814W 

F702W 

F555W,F702W 

F555W,F702W 

F555W,F702W 

F555W,F702W 

F555W,F702W,F814W 

F218W,F336W,F791W,F606W 

F702W,F791W 

F555W 



600,280 
600,280 
300 
280 

600,280 

460,280,460 

1700 

600,300 

600,280 

600,280 

600,300 

460,280,460 

2400,900,440,500 

280,750 

600 



Note. — Optical (WFPC2) images for each sample were obtained from the MAST (or 
from [http://acs.p ha. jhu.edu/ ^martel / ). The table indicates the proposal ID for each im- 
age used, along with the filter for the observation, and the integration time. [The complete 
version of this table is in the electronic edition of the Journal. The printed edition contains 
only a sample.] 



Table 7. Companion Sources 





Source 


JV 


C^.J2000 


^.72000 


D/" 


H 


PWHM 


Elong. 


Ellip. 


Sh 


TT 


M 


PM 


Maj 


Min 


Pt 


.Tot 


DD 


3C 15 


3 


00 37 04.05 


-01 09 05.5 


16.57 


19.21 


8.60 


1.26 


0.21 












Y 






Y 


3C 17 


3 


00 38 20.75 


-02 07 41.8 


6.73 


22.43 


7.15 


2.55 


0.61 








Y 




Y 








3C 20 


> 9 


00 43 09.67 


-1-52 03 30.5 


18.31 


17.80 


2.20 


1.04 


0.04 












Y 


Y 






3C 28 


> 9 


00 55 50.24 


-1-26 24 45.5 


22.65 


18.84 


4.31 


1.04 


0.04 












Y 




Y 




3C 29 


1 


00 57 35.33 


-01 23 37.4 


23.63 


24.04 


1.98 


1.47 


0.32 




















3C 31 


3 


01 07 24.44 


-1-32 24 44.0 


15.67 


21.20 


2.18 


1.28 


0.22 








Y 


Y 






Y 




3C 33 



















Y 












Y 




3C 33.1 


1 


oi 09 44.11 


+73 11 54.4 


4.49 


17.51 


4.87 


1.99 


0.50 








Y 


Y 










3C 35 


8 


01 12 01.93 


-1-49 28 44.2 


15.93 


17.81 


2.08 


1.06 


0.06 








Y 




Y 


Y 






3C 52 


> 9 


01 48 28.08 


-1-53 32 44.8 


6.60 


18.66 


5.12 


1.40 


0.29 












Y 


Y 


Y 




3C 61.1 


> 9 


02 22 36.07 


-1-86 19 07.7 


11.73 


18.10 


5.47 


1.12 


0.10 




Y 




Y 


Y 


Y 


Y 


Y 




3C 66B 


4 


02 23 11.54 


-H42 59 23.1 


16.55 


19.52 


2.60 


1.07 


0.06 








Y 


Y 




Y 




Y 


3C 71 


> 9 


02 42 41.38 


-00 00 49.4 


17.70 


18.06 


3.23 


1.13 


0.11 












Y 


Y 






3C 75N 


4 


02 57 43.63 


-1-06 02 31.5 


18.89 


18.74 


12.54 


2.69 


0.63 








Y 


Y 




Y 






3C 76.1 


5 


03 03 15.19 


-1-16 26 24.2 


9.65 


19.99 


17.36 


2.24 


0.55 




Y 








Y 


Y 


Y 





Note. — Number and properties of companion sources for each SCR target: SCR number; Number of companions on NIC2 image; Coordinates of brightest 
companion source; Distance (arcsec) of brightest companion from radio source; H-band absolute magnitude of brightest secondary; Diameter (FWHM) of 
secondary; Elongation (if resolved); Ellipticity (if resolved). The right-hand side of the table lists the presence of various types of artefact ('Y'): Shells; tidal 
tails; mergers; pre-mergers; major companion sources; minor companion sources; unresolved point-like companion sources; jets; dust disks. See the main test for 
a fuller description of each category. A question mark implies the presence of the artifact is suggested but cannot be confirmed categorically. [The complete 
version of this table is in the electronic edition of the Journal. The printed edition contains only a sample.] 
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3c66b 3clll 3c449 




0123 0123 0123 

r/arcsec r/arcsec r/arcsec 



Fig. 17. — 2DM radial profiles for the three test cases. Solid lines indicate best-fit 2DM 
model. Dashed lines indicate the best-fit nuclear component, with the dotted lines indicating 
the PSF scaled to match the nuclear flux of the source. 



A. 2. Synthetic galaxy tests 

We constructed simple synthetic galaxy models using Galfit (in non-optimization mode), 
and the parameters recovered from the three real test cases (3C 66B; 3C 111; 3C449). To 
these idealized galaxy models, varying amounts of noise were added, and different PSF's 
convolved with the data, in order to produce a synthetic NICMOS observation of an idealized 
galaxy. These synthetic observations were then modelled using 2DM. We began with the 
simplest case of the original models convolved with the PSF's, with no noise or background 
added. 2DM was then run just once to see how reliably it could recover the original parameter 
set used to define the synthetic observation. Results are presented for this "quick" run in 
table [TOl 2DM is found to recover the properties of the Galfit models to within ~ 10%. Note 
that in normal use, 2DM is run repeatedly to refine these fitting results. However, we adopt 
10% as a reasonable estimate for the systematic errors inherent in fitting non-ideal galaxies 
with idealized forms. 

Next we added in an additional random "readnoise" component to the simulated galax- 
ie s, and reran t he G alfit minimization, starting from a generic starting point, as described 



m 



Floyd et al.l (120041 ). Recovered parameters are presented in table [TTl The parameters are 



recovered almost perfectly from each galaxy with read-noise of 1 — lOe , plus Poisson. 



A. 3. Working Conclusions 

The results from 2DM and Galfit agree typically to within ~ 10% on morphological 
parameters, ~ 2% on host flux. However, there is a degree of variation, depending in 
particular on the treatment of errors in the central regions. This is in part due to the different 
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convolution algorithms applied. For highly luminous objects this leads to a disagreement in 
some host properties, while for all sources, there is still a large margin of error in the nuclear 
luminosity. 
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Table 8. Host galaxy R — H colors 



Source 




H 


R - H 


3cl5 


16.84 


13.55 


3.29 


3c29 


16.53 


12.89 


3.64 


3c31 


13.78 


11.18 


2.60 


3c35 


16.84 


13.68 


3.16 


3c66B 


14.38 


11.72 


2.66 


3c75N 


14.22 


12.07 


2.15 


3c83.1 


13.92 


11.53 


2.39 


3c84 


13.50 


12.45 


1.05 


3c88 


15.92 


12.23 


3.68 


3c98 


15.26 


13.03 


2.23 


3clll 


17.12 


14.83 


2.29 


3cl98 


17.64 


15.58 


2.06 


3c227 


16.21 


15.61 


0.60 


3c236 


17.62 


14.54 


3.08 


3c264 


14.17 


12.55 


1.62 



Note. — Approximate R — H colors 
for the sources based on ii-band fluxes 
at a 20 kpc aperture. 



Ref erences. — iMartel et al.| 
il999l) . [The complete version of this 
table is in the electronic edition of the 
Journal. The printed edition contains 
only a sample.] 
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Table 9. Galfit/2DM results: Real test-case galaxies 



Code 


Mag(liost) 


Re/kpc n{= 




b/a 


e 


Mag(nuc) 




3C 66b 


Galfit 


12.47 


2.93 


1.64 


0.92 


-79.4 


19.40 


2DM(nosamp) 


12.52 


2.67 


1.55 


0.91 


-86.5 


19.29 


2DM(samp) 


12.47 


2.81 


1.67 


0.91 


-85.9 


18.64 


3C 111 


Galfit 


14.40 


2.26 


6.98 


0.72 


16.8 


15.05 


2DM(nosamp) 


14.58 


2.24 


4.28 


0.63 


17.4 


15.04 


2DM(samp) 


14.71 


2.95 


3.75 


0.63 


17.2 


14.71 


3C 449 


Galfit 


12.48 


1.90 


1.36 


0.84 


-6.4 


21.34 


2DM(nosamp) 


12.51 


1.81 


1.30 


0.86 


-7.7 


20.46 


2DM(samp) 


12.49 


1.84 


1.37 


0.85 


-7.2 


20.52 



Note. — Comparison of results of fitting the three test case objects with 
Galfit and 2DM. Position angles are in degrees anti-clockwise from North. All 
magnitudes are uncorrected AB mag's from the PHOTFNU header keyword. 
Scale-lengths are half-light radii on the semi-major axis. The centre of the 
galaxy and nucleus are fixed at the same position in galfit as in 2DM, in order 
to force a comparison based on the same number of parameters. For 3C 111 
this posed problems (see text) so we present the galfit models both with and 
without a fixed centre. 
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Table 10. Galfit/2DM results: Synthetic Galaxies 



Code 


Mag(host) 


Re/kpc n(= 1/13) h/a 


e 


Mag(nuc) 




A (Simulated 3C 66B) 


Galfit 
2DM 


12.47 
12.38 


2.93 1.64 0.92 
3.15 1.89 0.90 


-79.4 
-72.4 


19.40 
19.87 


B (Simulated 3C 111) 


Galfit 
2DM 


14.40 
14.56 


2.26 6.98 0.72 
2.15 3.64 0.74 


16.8 
13.5 


15.05 
15.34 


C (Simulated 3C 449) 


Galfit 
2DM 


12.48 
12.49 


1.90 1.36 0.84 
1.87 1.34 0.84 


-6.4 
-6.3 


21.34 
21.22 



Note. — Comparison of results of fitting the three synthetic galax- 
ies, constructed with galfit using the parameters deduced for the three 
real test cases. Galfit recovers the parameter values to within 10% in 
a single run (no added noise). 
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Table 11. Response of Galfit to synthetic galaxies plus noise 



Galaxy 


Mag(host) Re/kpc n{= 




b/a 


e 


Mag(nuc) 




A (Simulated 3C 66b) 


Original 


12.47 2.93 


1.64 


0.92 


-79.4 


19.40 


Recovered(le~Poisson) 


12.47 2.93 


1.64 


0.92 


-79.4 


19.40 


Recovered (10e~Poisson) 


12.47 2.93 


1.64 


0.92 


-79.4 


19.41 


B (Simulated 3C 111) 


Original 


14.40 2.26 


6.98 


0.72 


16.8 


15.05 


Recovered (le~+Poisson) 


14.40 2.26 


6.93 


0.72 


16.8 


15.05 


Recovered (10e~+Poisson) 


14.39 2.26 


7.01 


0.72 


16.8 


15.06 


C (Simulated 3C 449) 


Original 


12.48 1.90 


1.36 


0.84 


-6.4 


21.34 


Recovered (le~+Poisson) 


12.48 1.89 


1.36 


0.84 


-6.4 


21.32 


Recovered (10e~+Poisson) 


12.48 1.89 


1.36 


0.84 


-6.4 


21.33 



Note. — Galfit results for fits to three synthetic test case galaxies constructed with 
varying degrees of readnoise (1 and 10 e~) plus Poisson noise. The three test cases are 
based on the Galfit models for 3C 66B, 3C 111, and 3C 449. The top line of each row 
gives the true parameters for the synthetic galaxy. The subsequent rows show the results 
recovered by running galfit. Each run was initiated at a generic starting point, not tailored 
to the individual cases. 



